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Abstract 

We present constraints on violations of Lorentz Invariance based on Lunar Laser Ranging (LLR) 
data. LLR measures the Earth- Moon separation by timing the round-trip travel of light between the 
two bodies, and is currently accurate to a few centimeters (parts in 10 11 of the total distance). By 
analyzing archival LLR data under the Standard-Model Extension (SME) framework, we derived 
six observational constraints on dimensionless SME parameters that describe potential Lorentz- 
violation. We found no evidence for Lorentz violation at the 10 -6 to 10 -11 level in these parameters. 

PACS numbers: 06.30. Gv, 11.30.Cp, 11.30.Qc, 95.55.Pe, 96.20.-n 
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Lorentz symmetry, the idea that physical laws take the same form in any inertial frame, ir- 
respective of the orientation or velocity of the frame, underpins the Standard Model and both 
Special and General Relativity. Attempts to quantize gravity have resulted in theories that 
allow for violations of Lorentz symmetry l|. For example, spontaneous Lorentz-symmetry 
breaking is possible in certain string theories At the Planck Energy, ^JTuFJg ~ 10 19 GeV, 
the Standard Model and General Relativity (GR) have comparable influence j^. Lorentz 
symmetry may not hold in that regime. Although no existing experiment can probe these 
energies, a Lorentz symmetry violation may also manifest itself at much lower energies 
where existing experiments do have sensitivity 4J. At present, no experiment has detected 
a violation of Lorentz symmetry. 

The Standard-Model Extension (SME) is a generalized effective field theory that adds 
a complete set of Lorentz- violating terms to the minimal standard model 141. The theory 

n n 

of the gravitational sector of the SME was developed in |5j and |6J. The SME provides 
a framework to analyze and compare the results of Lorentz symmetry experiments just as 
the Parametrized Post Newtonian (PPN) framework allows for the direct comparison of 
tests of gravity [7]. Recent calculations of the observable consequences of Lorentz symmetry 
violations in the pure-gravity sector of the minimal SME 6| showed that existing Lunar Laser 
Ranging (LLR) data would be sensitive to a subset of the pure-gravity SME parameters. 
LLR measures the time of flight of photons between a ranging station on the Earth and 
retro-reflectors on the lunar surface 8|. This experiment has been ongoing since the Apollo 
astronauts landed on the lunar surface in 1969. Over the past 38 years the measurement 
precision has improved by more than two orders of magnitude and currently LLR data can 
be used to determine the orbit of the Moon around the Earth to a few millimeters, or a 
few parts in 10 » of the tota. range Q. In this Letter, we present constraints on six 
linearly-independent combinations of pure-gravity sector SME parameters using archival, 
centimeter-precision LLR data. 



;he lunar orbit 
6J for our LLR 



A Lorentz violation would manifest itself as oscillatory perturbations to 
In this section, we describe the convention which we have adopted from 
data analysis. To explore these perturbations, it is convenient to work in a Sun-centered 
celestial-equatorial frame, (see Figured]). In this system, the Earth's equator defines the 
XY-plsme and the Earth's spin angular momentum vector points along the +Z direction. 
The Earth orbits the Sun in a plane (the ecliptic) which is inclined to the XY-plane by an 
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FIG. 1: The Sun-centered celestial-equatorial plane. The equator of the Earth defines the XY- 
plane and the Earth's rotation axis is aligned with the Z-axis. The orbital plane of the Earth 
(the ecliptic) is inclined to the XY-plane by the obliquity angle, rj ~ 23.5°. The Earth is shown 
at a descending node (a Vernal Equniox) which also defines the origin of our time coordinate, t. 
Reprinted figure with permission from 6J. Copyright 2006 by the American Physical Society. 



angle t] (approximately 23.5 degrees). The Earth's orbit intersects the Xl^-plane at two 
points, called the ascending and descending nodes. In Figure [U the Earth is pictured at 
the descending node, on the negative X-axis (note the arrow that indicates the direction of 
the Earth's motion along its orbit). At this time, commonly called the Vernal Equinox, the 
Sun, from the point of view of the Earth, ascends through the Earth's equatorial plane. The 
origin of t, our time coordinate (called T in jg]), coincides with the presence of the Earth at 
the descending node. 

Figure [2] depicts the orbit of the Moon about the Earth. The longitude of the ascending 
node and the inclination of the lunar orbit are labeled a and j3, respectively, and r is 
the mean orbital radius of the Moon. The lunar orbital phase, 9, corresponds to a unique 
reference frame for an arbitrary choice of the t — Vernal Equinox. The eccentricity and 
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FIG. 2: Lunar orbital parameters: Here the Earth is shown translated to the center of the Sun- 
centered coordinate system. The lunar orbit is described by ro, the mean distance between the 
Earth and Moon, e (not labeled) the eccentricity of the lunar orbit, a, the longitude of the ascending 
node, (3, the inclination of the lunar orbital plane with respect to the Earth's equatorial plane and 
9, the angle, along the Lunar orbit, subtended by the ascending node line and the position of 
the Moon at t = 0. Reprinted figure with permission from [g]. Copyright 2006 by the American 
Physical Society. 

longitude of perigee are not specified in this description. 

In this coordinate system, Lorentz-violating perturbations to the Earth- Moon separation 
can be expressed in the SME framework using the following fourier series: 

$r S ME(t) = Yj \ A n cos (u n t + n ) + B n sin (u n t + 0„)] . (1) 
n 

The dominant contributions to SrsMEif) occur at the following four frequencies: u, 2cu, 
2uj — loq and Here (as in 6j) u is the mean lunar orbital (sidereal) frequency, c^o 

is the anomalistic lunar orbital frequency (perigee to perigee) and f2® is the mean Earth 
orbital (sidereal) frequency. The corresponding amplitudes (A n and B n ) for these frequencies 
are listed in Table [H The dominant perturbations to the lunar orbit are controlled by a 
set of six linearly-independent combinations of Lorentz-violating SME parameters, sllr = 
{(s 11 — s 22 ) , s 12 , s 02 , s 01 , sn e ,c, ^n e ,s} (see jfj] for a description of these parameters). 

At present the only published constraints on the gravitational sector SME parameters 
come from measurements of the perihelion shifts of Mercury and the Earth and from an 
argument based on the current alignment of the solar spin axis and the angular momentum 
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TABLE I: The leading-order amplitudes from Equation Q] and their associated phases (from fjj). 
M and 5m are the sum and difference of the Earth and Moon masses, while V© and vq are the 
mean Earth and Moon orbital velocities, respectively, normalized to the speed of light. The unitless 
parameters b\ and b 2 are functions of 77, a, (3, ro, S~2© and the Earth's quadrupole moment. 

SME Parameter Partial Derivatives 

3( 5 n-^) = -f§ cos (2ut + 20) - j^l^ cos [(2a ; - W(J ) t + 20] 

5,5r/9s 12 = -f sin (2wt + 20) - 8( ^ e _ r ° o) sin [(2w - u> ) t + 29} 

d5r/d^ = -^^co S (ut + e) 

d5r/d^ = ^^ S m(u;t + 9) 

d5r/ds Q(BC = V®r (g.) cos (fi e <) 

<9<5r/asn es = ^©r (|y sin (fi e t) 
TABLE II: The partial derivatives of the SME perturbation to the lunar range with respect to 
each sllr parameter. 

vector of the planetary orbits jfj. This paper provides the first LLR constraints on SME 
parameters. 

The fundamental LLR observable is the travel time of light (usually a pulsed laser beam) 
that is propagated from a transmit station on the Earth to retroreflector arrays on the Lunar 
surface and is reflected back to a receive station on the Earth. Typically the receive station 
is the same as the transmit station. The return signal from the Moon is weak - typically 1 
to 5 photons per minute for the data used here - and single-photon counters are requred to 
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detect and time-tag lunar reflection events. 

LLR data have historically been presented as "normal points" which combine a series of 
single-photon lunar reflection events to achieve a higher signal-to-noise ratio measurement 
of the lunar range at some characteristic time for that data series. In this analysis, we 
used 14,401 normal points spanning September 1969 through December 2003, taken from 
the public LLR data archive [ill ]. The bulk of these normal points were generated by 



the McDonald Laser Ranging Station in Texas, USA |X2 1 and the Observatoire de la Cote 

n 

d'Azur station in Grasse, France [131 ] . A single normal point is typically generated from 
an observation spanning 5 to 20 minutes during which anywhere from a few to a hundred 
photons are collected. The lunar ranges reported in these normal points incorporate station- 
specific hardware corrections. 

In order to analyze LLR data, one relies on detailed models of the solar system ephemeris 
and dynamics. To our knowledge, there are only a few such models in existence capable of 



this analysis 



is fl, Q, 



151 ] . None of these codes explicitly include the SME framework in their 



equations of motion for solar system bodies. Nevertheless it is still possible to constrain the 
sllr parameters. 

We used the Planetary Ephemeris Program (PEP) [ijj] to extract SME parameter con- 
straints from LLR observations. PEP uses its ephemeris and dynamics model, along with 
an inital set of model parameters, to compute a range prediction and the partial derivatives 
of range with respect to each model parameter at the time of each normal point. It then 
performs a weighted, linear least-squares analysis to calculate adjustments to the model pa- 
rameters in order to minimize the difference between the observations and the model. Were 
the lunar range model linear in all parameters, this method would be exact. 

If one is concerned about non-linearities, one can solve for model parameters and then 
re-integrate the equations of motion, iterating until the parameter estimates converge. Over 
the past several decades the traditional (i.e. non-SME) analyses have done just that, re- 
sulting in agreement between model and data at the few centimeter level. Current model 
parameter values are therefore highly refined, and the weighted least-squares analysis sits 
firmly in the linear regime. As a result, it is not necessary to iterate when estimating new 
model parameters. Because the lunar range model is linear in the sllr parameters (see 
Equation [1] and Tabled]), the inclusion of these parameters in the analysis preserves linear- 
ity, as confirmed by the small adjustments to non-SME parameters seen in our solutions. 
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Performing an iterative solution for SME parameters requires the inclusion of the SME terms 
in the equations of motion. 

We computed the partial derivatives of lunar range with respect to each sllr parameter 
and provided this information to PEP prior to solving for the best-fit parameter adjustments. 
This approach is equivalent to explicitly including the sllr parameters in the equations of 
motion and setting their a priori values to zero (in which case SrsME{t) = so there is no 
SME contribution to the lunar orbit). We therefore treated any Lorentz violation as a small 
perturbation to a known orbit. The terms in the covariance matrix quantify any "crosstalk" 
between SME parameters and other fitted quantities. 

The solar system is complex. When modeling the expected light travel time between an 
LLR station on the Earth and a reflector on the lunar surface, one must account not only for 
the gravitational perturbations from the eight planets and Pluto, but also those of planetary 
satellites, asteroids, asphericities in the Sun, Earth and Moon, as well as various relativistic 
and non-gravitational effects (for a more complete description of relevant physical effects, 
see As a result there are many hundreds of parameters that have influence on the 

Earth-Moon range time. Not surprisingly, there are parameter degeneracies and LLR data 
alone cannot determine all of these parameters. Therefore LLR-only analyses suffer from 
systematic uncertainties in model parameter estimates and the formal error reported by the 
least-squares solution will underestimate the true model parameter uncertainties. These 
systematic uncertainties can dominate the formal errors. If auxiliary solar system data is 
included in the analysis (e.g. planetary radar ranging) then the number of model parameters 
grows and new parameter degeneracies arise. Having chosen to perform our analysis using 
LLR data alone, we accounted for the underestimation of parameter uncertainties by inflating 
the formal parameter uncertainties from our least-squares solution by a uniform factor which 
we label F . This is numerically equivalent to a uniform scaling of the uncertainty associated 
with each normal point by F. 

In order to determine the F factor, we performed an analysis of the LLR data in which 
we froze all sllr parameters at zero, but allowed the PPN parameters, f3 and 7 to vary. We 
know from other experiments 0, [if], Q, 3 that (3 and 7 are consistent with their GR value 
of unity to within a part in 10 3 or better. We found that we must scale our uncertainties on 
j3 and 7 by F = 20 to be in accord with these earlier results. 

We then estimated the values of the set of sllr parameters while holding the PPN 
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and the values derived 



TABLE III: The predicted sensitivity to each sllr parameter (from 
in this work including the realistic (scaled) uncertainties (Fa) with F = 20. In this analysis, the 
PPN parameters were fixed at their GR values. The SME parameters are all within 1.5 standard 
deviations of zero. We see no evidence for Lorentz violation under the SME framework. 

parameters at their GR values (f3 — 7 = 1). The resulting parameter values and their 
realistic errors (the formal errors scaled by the F = 20 factor) are reported in Table IIHI All 
SME parameters are within one and a half standard deviations of zero and are within an 
order of magnitude of the sensitivity expected from 1 cm precision LLR data as estimated in 
6|. There is no evidence for Lorentz violation in the lunar orbit. The model-data agreement, 
binned annually, is shown in Figure O 

To ensure that our analysis is sensitive to the signature of a Lorentz violation in the lunar 
range, we generated a perturbed LLR dataset with a hand-inserted Lorentz-violation signal 
corresponding to a value of s 11 — s 22 = 9 x 10~ 10 (a factor of 10 larger than the uncertainty for 
this parameter). Our analysis recovers this signal. We found s 11 — s 22 = (10 ± 0.9) x 10" 10 
with the other sllr parameters unchanged. 

In conclusion, we have analyzed over 34 years of archival LLR data using PEP to derive 
constraints (see Table IIIII) on six linearly-independent Lorentz-violation parameters in the 
pure-gravity sector of the SME. These are the first LLR-based SME constraints. Using the 
same LLR dataset, these constraints could be improved by performing a simultaneous fit 
with auxiliary solar system data to help break parameter degeneracies and thereby reduce 
the systematic error budget. Based on the solutions for f3 and 7 obtained in this work, we 
estimate that the SME parameter constraints could be improved by a factor of 5-10 with a 
joint fit. In addition, we are involved in the Apache Point Observatory Lunar Laser-ranging 
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FIG. 3: The annual RMS residual between the LLR data and our best-fit model for the lunar 
range. The residual RMS in 1969 is over 300 cm. We omitted this point from the plot for clarity, 
but the two data points from that year were included in the analysis. Over this timespan, the 
potential Lorentz-violating signals would have undergone at least 34 cycles. 

Operation project (APOLLO), a next-generation LLR station that is currently collecting 



millimeter-precision lunar range data 



101 ] . We plan to improve the SME parameter con- 



straints by incorporating APOLLO data into our analysis. 

This project developed out of conversations with several colleagues including E. Adel- 
berger, Q. Bailey, J. Davis, A. Kostelecky, J. Moran, T. Murphy, I. Shapiro and M. Zaldar- 
riaga. We would also like to thank the LLR observers and the ILRS for the data collection 
and normal point distribution. J. B. R. B. acknowledges financial support from the ASEE 
NDSEGF, the NSF GRFP and Harvard University. Generous financial support was provided 
by the National Science Foundation (Grant No. PHY-0602507). 
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